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Abstract 

It is shown that, if the "new neutrino" imphed by the Reactor Neutrino 
Anomaly exists and is in fact characterized by the suggested relatively high 
mass squared difference and reasonably large mixing angle, it should clearly 
reveal itself in the oscillometry measurements. For a judicious neutrino source 
the "new oscillation length" L14 is expected shorter than 3m. Thus the 
needed measurements can be implemented with a gaseous spherical TPC of 
modest dimensions with a very good energy and position resolution, detecting 
nuclear recoils following the coherent neutrino-nucleus elastic scattering. The 
best candidates for oscillometry, yielding both monochromatic neutrinos as 
well as antineutrinos, are discussed. A sensitivity in the mixing angle ^14, 
sin^ (26'i4)=0.1 (99%), can be reached after a few months of data handling. 

Keywords: sterile neutrinos, oscillometry,neutral currents, spherical 
gaseous TPC. 



1. Introduction. 

A recent analysis of the Reactor Neutrino Anomaly (RNA) jl[ led to a 
challenging claim that this anomaly can be explained in terms of a new fourth 
neutrino with a mass squared difference much larger than one encounters in 
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neutrino oscillations. In fact assuming that the neutrino mass eigenstates 
are non degenerate one finds 

a2 i2 21 a2 i2 2i /\2 /\ 

Am24 = — m^l f» Am;^4 = — > 1.5(eV) . (1) 
and a mixing angle 

sin^ 2^14 = 0.17 ±0.1(95%). (2) 

It is obvious that this new neutrino should contribute to the oscillation 
phenomenon. In the present paper we will assume that the new netrino is 
sterile, that is it does not participate in weak interaction. Even then, how- 
ever, it has an effect on neutrino oscillations since it will tend to decrease 
the elecron neutrino flux. This makes the analysis of oscillation experiments 
more sophisticated. In all the previous experiments the oscillation length is 
much larger than the size of the detector. So one is able to see the effect only 
if the detector is placed in the right distance from the source. It is, however, 
possible to design an experiment with an oscillation length of the order of 
the size of the detector, as it was proposed in [2[,j3[. This is equivalent to 
many standard experiments done simultaneously. 

In a previous paper [3] we have studied the possibility of investigating the 
low energy neutrino oscillations by the measurement of the electron recoils 
following neutrino electron scattering. In the present study we will explore 
the neutral current interaction to measure nuclear recoils. The main require- 
ments are the same as given previously ^ with some modifications pertinent 
to the neutrino nucleus interaction. More specifically: The neutrinos should 
have low energy so that the oscillation length is smaller than the size of the 
detector. At the same time it should sufficiently high so that the neutrino- 
nucleus elastic scattering yields recoils above threshold with a sizable cross 
section. A monoenergetic neutrino source is preferred, since it has the ad- 
vantage that some of the features of the oscillation pattern are not affected 
as they might be by the averaging over a continuous neutrino spectrum. 
Antineutrino sources with a relatively high energy can also be employed. 
The "wave form" of the oscillation is not much different than that of the 
monochromatic source, since only the small portion of the spectrum at the 
high energy end becomes relevant. The lifetime of the source should be suit- 
able for the experiment to be performed. Clearly a compromise has to be 
made in the selection of the source. 

In this article we will show that, unlike the standard neutrino case, for a 
sterile neutrino one can observe neutrino diasapperance oscillations via the 
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neutral current interaction. Furthermore the aim of this article is to show 
that the existence of a new fourth neutrino can be verified experimentally 
by the direct measurements of the oscillation curves for the monoenergetic 
neutrino-nucleus elastic scattering. It can be done point-by-point within the 
dimensions of the detector, thus providing what we call neutrino oscillometry 



2. Neutrino oscillations and neutral current detection 

Suppose in addition to the three standard neutrinos we have a 4th sterile 
neutrino. In the neutral current detection all contributing neutrinos have the 
same cross section, say a. Let us suppose that we initially have electronic 
neutrinos Ve- Then we distinguish the following cases: 

• All four neutrinos are active. 
Then 

^tot = (^(^e Ve) + P{ye ^ + Pi^e ^ ^r) + P{l^e ^ ^i)) O , 

(3) 

but 

P(l/e ^ I/e) = 1 - {P{V^ U^) + P(z/e ^ U^) + P(z/e ^ U^)) , (4) 

i.e. 

^tot = (5) 

no oscillation is observed. 

• The fourth neutrino is sterile. 
Then 

^tot = (Pi^e ^ I'e) + P{Ve ^i?) + P{Ve ^r)) C^, (6) 

i.e. the sterile neutrino does not contribute. Eq. HJ however, is still 
valid (neutinos are lost from the flux). Thus 

^tot = (l-^('^e^^^4))cT. (7) 

If, in addition, the new oscillation length is much smaller than the other 
two, one finds: 

^tot ^ (^■'" ~ ^^^^ ^^^^ ^' 
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It is worth comparing the neutral current situation with neutrino-electron 
elastic scattering previously considered Then for the standard neutrinos 
we have 

^tot = Pi^e ye)(y + {Pi^e ^ + Pi^e ^ ^r)) o', (9) 

where o is the (i/g, e) cross section, while a' is the cross section for the other 
two flavors, (i/q,, e), a = fi,T,. Furthermore 

P(z/e ^ Z/e) = 1 - (P(z/e ^ U^) + P{Ue ^ ^^r)) • (10) 

If, in addition, the new oscillation length is much smaller than the other two, 
one finds: 

atot ~ (l - xiE,) sin' 2^13 sin' 7r£-^ a, , x{K) = (^1 - , (11) 
i.e. one observes oscillations, since j2l,|3| o' ^ a. 

Furthermore, if there exists an additional sterile neutrino, we find that the 
corresponding oscillation cross section in (z/, e) scattering is still given by Eq. 

3. The differential and total cross section 

The elastic neutrino-nucleus scattering due to the neutral current inter- 
action has beenpreviously considered for the detection of sky 0] and 
earth neutrinos [Tl] of appreciably higher neutrino energies. It has never been 
considered in the context of neutrino oscillations, since for standard neutri- 
nos, as we have seen, oscillations in such a channel are not expected. The 
differential cross section for a given neutrino energy Ey can be cast in the 
form [2I: 

(T^, E;) = ^Hp^ {NyA)FUTA, Ey), (12) 

with 

FUTa, Ey) = F'(g') (1 + (1 - - , (13) 

where N is the neutron number and F{q^) = F{T]^ + 2A'mNTA) is the nuclear 
form factor. The effect of the nuclear form factor depends on the target. In- 
tegrating the total cross section of Eq. [T2] from = Tfh to the maximum 
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allowed by the neutrino energy we obtain the total cross section. The thresh- 
old energy Tt^ depends on the detector. Since the favorable neutrino energies 
are of order 1 MeV, the results are not sensitive to the nuclear form factor. 
They crucially depend on the detector threshold, since the energy of the re- 
coiling nucleus is quite low. This retardation becomes more severe, in the 
presence of quenching. 

Furthermore for a real detector the expected nuclear recoil events are 
quenched, especially at low energies. The quenching factor for a given detec- 
tor is the ratio of the signal height for a recoil track to that of an electron 
signal with the same energy. We should not forget that the signal heights 
depend on the velocity and how the signals are extracted experimentally. 
The actual quenching factors must be determined experimentally for each 
target. In the case of Nal the quenching factor is 0.05, while for Ge and Si it 
is 0.2-0.3. Thus the measured recoil energy is typically reduced by a factor 
of about 3 for a Si-detector, when compared with the electron energy. For 
our purposes it is adequate, to multiply the energy scale by an recoil energy 
dependent quenching factor, QfaciTA) adequately described by the Lindhard 
theory jsf, [oj. More specifically in our estimate of Qfacl^^) we assumed a 
quenching factor appropriate for a gas target of ^He which fits well the data 
lol | in the energy region of 2 to 50 keV : 



Ta 
IkeV 



n ~ 0.620 , ra ~ 0.070. (14) 



The quenching factor as given by Eq.( [HI) is exhibited in Fig. [11(a) for the 
energy of interest to us. We will assume that it is the same for all noble 
gases of interest in the present work. Due to quenching the threshold energy 
is shifted upwards, from Tth to (see Fig. [11(b)) 

The minimum neutrino energy required as a function of threshold is pre- 
sented in Fig. [21 Furthermore, even if the neutrino energy is above this 
minimum and the detection is allowed with a coherent cross section (oc A^^), 
the heavier the target, the smaller the recoiling energy is and the more effec- 
tive the retardation due to the threshold becomes. 

Once the usual neutrino-nucleus cross sections aA^E^, 0) are known one 
can show that, due to the oscillation, they depend on the distance of the 
obsevation point from the source. One finds that the neutrino disappearance 
as seen in the neutrino-nucleus cross section can be cast in the form: 



aA{E,,L) = aA{E„0) 



1 - sin2 2^i4sin2 ( 3.72451L 



me 
E. 



(15) 
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(a) (b) 
Ta keV 

Figure 1: The quenching factor as a function of the recoil energy of interest 
in the present work (a). Due to quenching the threshold energy for nuclear 
recoils is shifted upwards from Tth to T^, e.g. from 0.10 to 0.18 keV (b). 




(a) (b) 

T„,i„ keV 

Figure 2: We show the minimum neutrino energy required as a function of 
threshold without quenching (a) and with quenching in (b). From top to 
bottom for the targets of ^^^Xe, '^"Ar, ^°Ne and ^He. The threshold value is 
very crucial, especially for heavy targets. 



The obtained results crucially depend on the enegy threshold and the 
quenching factor. Even though the threshold achieved by the gaseous spheri- 
cal time projection counter (STPC) is impressive, 0.1 keV, since the neutrino 
recoiling energy is extremely low, in particilar for heavy targets, most of the 
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neutrino sources do not pass the test. Some of those which pass the test are 
included in table [TJ From these sources we will consider for further analysis 
those, which can be employed at present, namely ^^Ar, ^^Cr, ^^Zn and ^^P. 

Table 1: Proposed candidates for a new neutrino oscillometry at the spherical 
gaseous TPC. Tabulated nuclear data have been taken from other data 
have been calculated in this work (see the text for details. The mass of 
the source was assumed to be 0.1 Kg, except for ^^Ar and the antineutrino 
source, ^^P, for which we took 0.0043 kgr each). We note that the oscillation 
lengths for the new neutrino [ij , 12 1 estimated in this work are much smaller 



compared to those of all the earlier experiments. 



Nuclide 


Tl/2 




Li4 


K 




(d) 


(keV) 


(m) 


(s-i) 


37 Ar 


35 


811 


1.4 


1.8 X 


51 Cr 


27.7 


747 


1.2 


4.1 X 10^^ 


«5Zn 


244 


1343 


2.2 


3.0 X 


59Ni 


2.8 X 10'^ 


1065 


1.8 


1.1 X 


ii^Sn 


116 


617 


1.0 


3.7 X 


32 p 


14.3 


continuum 


^ 2.5 


5.0 X 10^6 



4. Results with monochromatic sources 

We will begin with the cross sections and proceed to the event rates. 

4.1. Coherent cross sections 

For a light target like He, a variety of sources pass the test, but one 
loses the benefit of large coherence. In comparing these results with those of 
neutrino-electron scattering we should keep in mind that the present cross 
sections explicitly contain the A^^ term due to coherence, in the case of 
neutrino-electron scattering the enhancement Z is not included in the cross 
section, but it has been contained in the electron densiity, which is Z times 
the number of nuclei per unit volume. 

For the sources and the targets that pass the test we present our results 
in FigsEB 
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(a) 



L — )■ meters 



(b) 



Figure 3: In panel (a) we show the total cross section in units of 
{G pTneY / {2tt) = 2.29 x 10~'*^m^ as a function of the distance of the re- 
coiling nucleus from the source in meters assuming a threshold of 0.1 keV. 
In panel (b) we show the same quantity in the presence of quenching. 
In both panels the solid, dotted and dotted-dashed curves correspond to 
sin^26'i4 = 0.27,0.17,0.07 respectively. The results were obtained for a Ne 
target with the ^^Zn source. 



4-2. The event rate 

We will consider a spherical detector with the source at the origin and 
will assume that the volume of the source is much smaller than the volume 
of the detector. The event rate dl between L and L + dL is given by: 

dl = NyHA— — :p^a{L,x) = N^nAdLa{L, x), (16) 
where N,y is the neutrino intensity and ua is the number of the target nuclei: 

where P is its pressure and Tq its temperature. Thus 

dl , , , s 

— = N,nAa{L,x) (18) 
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(a) (b) 
L — )■ meters 



Figure 4: In panel (a) we show the total cross section in units of 
{G FTUeY / {2tt) = 2.29 X 10~'*^m^ as a function of the distance of the re- 
coiling nucleus from the source in meters assuming a threshold of 0.1 keV. 
The results were obtained for a He target with the ^^Ar source. In panel 
(b) we show the same quantity for a ^^Cr source. In both panels the ef- 
fect of quenching was included. The solid, dotted and dotted-dashed curves 
correspond to sin^ 26*14 = 0.27, 0.17, 0.07 respectively. 



R^— = AHL,x), (19) 

where 

A = ^^RoKriA (20) 

Rq the radius of the target and a{L,x) is the neutrino-nucleus cross section 
in units of [G pTrieY / {^n) = 2.29 x 10~^^m^. The total number of events per 
unit length after running for time tr will be given by 

^ = A^a{L, x) (1 - e-*^/^ , A^ = ^riV.n^ (21) 

where r is the lifetime of the source. Integrating Eq. (I2T1) over L from to 
Rq we obtain the total number of events, which can be cast in the form: 

N = A + Bsm^2eu. (22) 
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The parameters A and B for some cases of interest are presented in table [2J 
We have included here the relevant results even for the pairs source-target 
not discussed further in the present work in the context of oscillometrjQ, just 
to give an estimate on the uncertainties expected in the extraction of sin^ 2^14 
from the total number of events 

Table 2: The parameters A and B entering the total event rate (for their 
definition see the text, Eq. fl22]) ). The source ^"^Bi is included to show the 
importance of the size of the energy of a monochromatic source, even though 
it is, at present, unrealistic (in the sense of production) to employ it. 



target- 


A (no quenching) 


B (no quenching) 


A (quenching) 


B (quenching) 


source 










■^^Ar- ^^P 


2.4 X 10^ 


-1.2 X 10^ 






40 Ar- 205Bi 


1.4 X 10^ 


-6.6 X 10^ 


4.2 X 10^ 


-1.8 X 10^ 


20Ne_32p 


8.8 X 10^ 


-4.6 X 10^ 


1.0 X 10^ 


-5.4 X 10 


20Ne- 65Zn 


2.9 X 10^ 


-1.6 X 10^ 


5.3 X 10^ 


-2.8 X 10^ 


20Ne_ 205Bj 


7.2 X 10^ 


-3.3 X 10^ 


3.8 X 10^ 


-1.7 X 10^ 


^He- 37Ar 


7.8 X 10 


-3.9 X 10 


3.6 X 10 


-1.8 X 10 


4He-"Cr 


8.7 X 10^ 


-4.1 X 10^ 


3.1 X 10^ 


-1.5 X 10^ 


^He- ^^Zn 


4.0 X 10^ 


-2.1 X 10^ 


3.3 X 10^ 


-1.8 X 10^ 


4He-205Bi 


4.6 X 10^ 


-2.0 X 10^ 


4.3 X 10^ 


-1.9 X 10^ 



The goal of the experiment is to scan the monoenergetic neutrino nucleus 
elastic scattering events by measuring the nuclear recoils as a function of 
distance from the neutrino source prepared in advance at the reactor/s. This 
scan means point-by-point determination of scattering events along the detec- 
tor dimensions within its position resolution. These events can be observed 
as a smooth curve, which reproduces the neutrino disappearance probability. 
It is worthwhile to note again that the oscillometry is suitable for monoener- 
getic neutrino, since it deals with a single oscillation length or Li4(see table 
[1]). This is obviously not for antineutrino, since, in this instance, one 



^ Due to lack of space the expected oscillation patterns are not provided in this paper. 
They can be obtained by communicating directly with the authors 
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extracts only an effective oscillation length. This could be a serious problem 
in neutrino electron scattering. In the case of nuclear recoils, as we will see 
below, with a judicious choice of the target-source pair, not much information 
is lost due to the folding, since only a narrow band in the high energy tail of 
the continuous neutrino energy spectrum contributes to nuclear recoils, even 
though the assumed threshold for noble gas targets is quite low, 0.1 keV. 

Table [T] clearly shows that the oscillation lengths for a new neutrino pro- 
posed in [l| are much smaller compared to those previously considered 
in connection with ^13. They can thus be directly measured within the di- 
mensions of detector of reasonable sizes. One of the very promising options 
could be the Spherical Time Projection Counter (STPC) proposed in M. If 
necessary, a spherical Micromegas based on the micro-Bulk technology |13| . 
which will be developed in the near future, can be employed in the STPC. A 
thin 50 micron polyamide foil will be used as bulk material to fabricate the 
detector structure. This detector provides an excellent energy resolution, can 
reach high gains at high gas pressure (up to 10 Atm) and has the advantage 
that its radioactivity level ij] should fulfill the requirements of the proposed 
experiment. 

In this spherical chamber with a modest radius, assumed to be 4m, the 
neutrino source can be situated in the center of the sphere and the detector 
for recoils is also placed around the source in the smaller sphere with radius 
r ^ 0.5 m. The chamber outside this small sphere is filled with a gas (a 
noble gas such as Ar, Ne, or He (preferably Ar), if the neutrino source is of 
sufficiently high energy). In the present work we assumed the gas is at room 
temperature under a pressure of 10 Atm. 

The nuclear recoils are guided by a Micromegas-detector [l5|,[l^. Such 
type of device has an advantage in precise position determination (better than 
0.1 m) and in detection of very low nuclear recoils in 47r-geometry (down to 
a few hundreds of eV, that well suits the nuclides from Tabled]). 

The results obtained in the presence of quenching are presented in Figs 



5. Antineutrino Sources 



The monochromatic neutrino sources considered above, unfortunately, 
have the disadvantage that the neutrino energy is lower than required to meet 
the experimental requirements for a neutral current (NC) detector with high 
neutron number. Thus one may have to resort to antineutrino sources paying 
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(a) (b) 
L — meters 

Figure 5: In panel (a) we show the number of events as a function of the 
distance of the recoihng nucleus from the source in meters assuming a thresh- 
old of 0.1 keV in the presence of quenching. The results were obtained for 
a He target with the ^^Ar source, i.e. Ai = 1.9m~^. In panel (b) we show 
the same quantity with the ^^Cr source, i.e. Ai = 28.8m~^. In both panels 
the target gas was enclosed in a sphere of 4m radius, at temperature T=300 
°K and under a pressure of 10 Atm. In both panels the solid, dotted and 
dotted-dashed curves correspond to sin^ 26*14 = 0.27, 0.17, 0.07 respectively. . 



the price of distortion of the form of the oscillation due to the integration 
over the antineutrino spectrum. The NC cross section for an ^*^Ar target with 
a threshold of 0.1 keV as a function of neutrino energy is shown in Fig. [71 
This cross section must be folded with the energy spectrum of the source. 
From an experimental point of view '^^P is the best antineutrino source (see 
tabled]). In principle, even though it has not been included in table [H ^°Sr 
can also be used as an antineutrino source. In the case of ^^P the normalized 
spectrum is exhibited in Fig. |8l One sees that, due to theshold effects only 
a portion of the spectrum can be exploited. Since the regions of large cross 
section have a small probability in the spectrum, the integrated cross section 
is not very large. In fact folding Eq. US] we obtain the results shown in Figs 
[UlfTUl From the thus obtained cross sections proceeding as above we obtain 
the differential number of events dN/dL exhibited in Fig. [11] Integrating 
over L the results of Fig. [TT] and the corresponding one without quenching, 
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(a) (b) 
L — meters 

Figure 6: The same as in Fig. Elfor a Ne target (a) and He target (b) with 
the ^^Zn source, i.e Ai = 22.6m~^. 




1.45 1.50 1.55 1.60 1.65 1.70 1.1 1.2 1.3 1.4 1.5 1.6 1.7 

(a) (b) 

Figure 7: The non oscillating part of the antineutrino cross section in units 
of {G FrrieY / {2tx) = 2.29 x 10~''^m^ as a function of the energy in MeV for a 
target ^°Ar (a) and ^^Ne (b), assuming a threshold of 0.1 keV. 



not shown here, we obtain the parameters A and B shown in table O 
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Ey ->MeV 



Figure 8: The normalized antineutrino spectrum following the beta decay of 
^^P. The vertical line indicates the space on its right allowed for ^°Ar (dotted 
line) and ^°Ne (dashed line) targets, assuming a threshold of 0.1 keV. 



6. Discussion 

We have seen that neutrino oscillometry will lead to a direct observation of 
the fourth sterile neutrino in electron neutrino disapperarnce experiments, if 
such a neutrino exists. The calculations and analysis shows that the gaseous 
STPC is a powerful tool for identification of a new neutrino in two ways 
i) by observing the electron recoils in neutrino-electron scattering , recently 
discussed by us 0], and ii) by measuring the nuclear recoils in the coherent 
neutrino-nucleus scattering via the neutral current interaction, discussed in 
this work. The latter method can utilize the wide experience obtained in 
other experiments attempting to measure nuclear recoils, like dark matter, 
(quenching factors, background minimization etc). 

We have seen that, since the expected mass-squared difference for this 
neutrino is rather high, the corresponding oscillation length is going to be 
sufficiently small for 1 MeV neutrino energy, so that it can be fitted into 
the dimensions of a spherical detector with the radius of a few meters. The 
neutrino oscillometry can be implemented in this detector with the use of 
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3 4 I 2 3 4 — ). 

(a) (b) 
meters 

Figure 9: The antineutrino cross section in units of (G'i?me)^/(27r) = 2.29 x 
10~^^m^, folded with the spectrum of Fig. [HI as a function of L in m for a 
target ^°Ar, assuming an energy threshold of 0.1 keV and ^^P source, without 
the quenching factor (a) and including quenching (b).The solid, dotted and 
dotted-dashed curves correspond to sin^ 26^14 = 0.27, 0.17, 0.07 respectively. 



the intense neutrino sources which can be placed at the origin of the sphere. 
The gaseous STPC with the Micromegas detection has a big advantage in 
the 47r-geometry and in very good position resolution (better than 0.1 m) 
with a very low energy threshold (~ 100 eV). 

The most promising candidates for oscillometry have been considered. 
These are the monochromatic sources '^^Ar, ^^Cr and ^^Zn as well as ^^P for an 
antineutrino source. As an example, for the ^^Zn source the sensitivity to the 
mixing angle 6'i4 is estimated as sin^ (26*14) = 0.1 with the 99% of confidence, 
arising from the total number of events collected during only a few months 
of data handling. The observation of the expected characteristic oscillometry 
curve will provide a much more precise information on the oscillation length 
and thus constitute a definite manifestation of the existence of a new type of 
neutrino as very recently proposed by the analysis of the reactor antineutrino 
anomaly. 

Unfortunately, however, one cannot fully benefit from the large size of 
the coherent cross section for large neutron number N of the target. The 
reason is that larger N implies larger nuclear mass number A, i.e. lower 
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(a) 



L — )■ meters 



(b) 



Figure 10: The same as in Fig. |9] for the target ^°Ne and the ^^P source , 
assuming a threshold of 0.1 keV, without the quenching factor (a) and in- 
cluding quenching (b). At the position of the second peak, we see a distortion 
of the "wave form", which, without quenching, is bigger due to the fact that 
the allowed neutrino energy range is, in this case, wider. 



recoil energy. Thus, even though the gaseous STPC detector has already 
achieved the impressive threshold of 0.1 keV, this is still high enough to 
reduce the cross sections for neutrinos in the MeV range for heavy targets or 
even exclude these targets. It is clear that the accuracy of the extraction of 
the above parameters will substantially increase only if the threshold can be 
further reduced, since it does not seam feasible to utilize a monochromatic 
neutrino source with energy higher than that of 1343 keV associated with 
^^Zn. 
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